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Abstract 

In this paper, we consider the development of a class of one-step hybrid third derivative 
block method with three off-grid points for the direct solution of initial value problems 
of second order Ordinary Differential Equations. We adopted method of interpolation and 
collocation of power series approximate solution to generate the continuous hybrid linear 
multistep method, which was evaluated at grid points to give a continuous block method. 
The discrete block method was recovered when the continuous block method was evaluated 
at selected grid points. The basic properties of the method was investigated and was found 
to be zero-stable, consistent and convergent. The efficiency of the method was tested on 
some stiff equations and was found to give better approximation than the existing method, 
which we compared our result with. 
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Introduction 

This paper solves second order initial value problems in the form 

y" = f(x,y(x),y\x)) , y(x 0 ) = y 0 , y'(x 0 ) = y' 0 (l) 

where f is continuous within the interval of integration. Solving higher order derivatives 
method by reducing them to a system of first-order approach involves more functions 
to evaluate which then leads to a computational burden as in [1-3], The setbacks of this 
approach had been reported by scholars, among them are Bun and Vasil'yer and Awoyemi 
et al [4,5], 

The method of collocation and interpolation of the power series approximation to 
generate continuous linear multistep method has been adopted by many scholars; among 
them are Fatunla, Awoyemi, Olabode, Vigor Aquilar and Ramos, Adeniran et al, Abdelrahim 
et al, Mohammad et al to mention a few [5-11], Block method generates independent 
solution at selected grid points without overlapping. It is less expensive in terms of number 
of function evaluation compared to predictor corrector method, moreover it possess the 
properties of Runge Kutta method for being self-starting and does not require starting 
values. Some of the authors that proposed block method are: [12-17], 

In this paper, we developed a on e-step hybrid third derivative method with three 
offgrid, which is implemented in block method. The method is self-starting and does not 
require starting values or predictors. The implementation of the method is cheaper than 
the predictor-corrector method. This method harnesses the properties of hybrid and third 
derivative, this makes it etcient for stin problems. 

The paper is organized as follows: In section 2, we discuss the methods and the materials 
for the development of the method. Section 3 considers analysis of the basis properties of 
the method which include convergence and stability region, numerical experiments where 
the etciency of the derived method is tested on some numerical examples and discussion of 
results. Lastly, we concluded in section 4. 

Derivation of the Method 

We consider a power series approximate solution of the form 
2 r+s- \ ( v _ v 

y( x )= E a i [2) 

y=0 V 

where r=2 and s=5 are the numbers of interpolation and collocation points respectively, 
is considered t o be a solution to (1). 

The second and third derivative of (2) gives 
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Substituting (3) into (1) gives 
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Collocating (4) at all points x„ + ,,r = 0^-jl and Interpolating 
Equation (2) at x j-S 0,— . gives a system of non linear equation 
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Solving (5) for a t ' S using Gaussian elimination method, gives a 
continuous hybrid linear multistep method of the form 


y{ x )= Z «sfai +A! 


I/?/, +/9/ t 
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Differentiating (6) once yields 
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Equation (6) is evaluated at the non-interpolating points and (7) 
at all points, produces the following general equations in block form 
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Equation (6) is evaluated at the non-interpolating points 
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Multiplying equation (8) by the inverse of (A) gives the 
hybrid block formula of the form 
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Analysis of Basic Properties of the Method 
Order of the Block 

According to [7] the order of the new method in Equation (8) is 
obtained by using the Taylor series and it is found that the developed 
method has a uniformly order eleven, with an error constants vector 
of: 

~6.1829x 1CT 14 ,1.752x1 (T 13 ,3.0454 x 1(T 13 , 4.8542x 10~ 13 ,T 

C = 

" [_ 4.1791xl0~ 13 ,4.8542x 10~ 13 ,5.5292x 10~ 13 ,9.7083x 10~ 13 _ 

Consistency 

The hybrid block method [4] is said to be consistent if it has an 
order more than or equal to one. 

Therefore, our method is consistent. 

Zero Stability of Our Method 

Definition: A block method is said to be zero-stable if as h —> 0 
, the root z j ,i = \(Y)k of the first characteristic polynomial 
p(z) = 0 that is 
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p(z) = z 6 (z-l) 2 =0, 

Hence, our method is zero-stable. 

Regions of Absolute Stability (RAS) 

Using Mat Lab package, we were able to plot the stability region 
of the block methods (see figure 1 below). 

The stability polynomial for K=1 with three offstep point using 
Scientific Workplace software package we have in the following one: 
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x-values 

Exact Solution 

Computed Solution 

Error in our method 

Error in [13] 

0.100 

0.90483741803595957316 

0.90483741803595957245 

7.100E(-19) 

1.054712E(-14) 

0.200 

0.81873075307798185867 

0.81873075307798185795 

7.200E(-19) 

1.776357E(-14) 

0.300 

0.74081822068171786607 

0.74081822068171786534 

7.300E(-19) 

2.342571E(-14) 

0.400 

0.67032004603563930074 

0.67032004603563930001 

7.300E(-19) 

2.797762E(-14) 

0.500 

0.60653065971263342360 

0.60653065971263342285 

7.500E(-19) 

3.130829E(-14) 

0.600 

0.54881163609402643263 

0.54881163609402643185 

7.800E(-19) 

3.397282E(-14) 

0.700 

0.49658530379140951470 

0.49658530379140951390 

8.000E(-19) 

3.563816E(-14) 

0.800 

0.44932896411722159143 

0.44932896411722159058 

8.500E(-19) 

3.674838E(-14) 

0.900 

0.40656965974059911188 

0.40656965974059911099 

8.900E(-19) 

3.730349E(-14) 

1.00 

0.36787944117144232160 

0.36787944117144232065 

9.500E(-19) 

3.741452E(-14) 


Table 1: Comparison of the proposed method 


x-values 

Exact Solution 

Computed Solution 

Error in our method 

Error in [16] 

.0.100 

-0.10517091807564762480 

-0.10517091807564762481 

1.000E(-20) 

8.326679(-17) 

0.200 

-0.22140275816016983390 

-0.22140275816016983391 

1.000E(-20) 

2.775557(-16) 

0.300 

-0.34985880757600310400 

-0.34985880757600310397 

-3.000E(-20) 

5.551115(-16) 

0.400 

-0.49182469764127031780 

-0.49182469764127031779 

-1.000E(-20) 

9.436896(-16) 

0.500 

-0.64872127070012814680 

-0.64872127070012814679 

-1.000E(-20) 

2.109424(-15) 

0.600 

-0.82211880039050897490 

-0.82211880039050897478 

-1.000E(-19) 

3.219647(-15) 

0.700 

-1.01375270747047652160 

-1.01375270747047652150 

-1.000E(-19) 

4.440892(-15) 

0.800 

-1.22554092849246760460 

-1.22554092849246760440 

-2.000E(-19) 

5.995204(-15) 

0.900 

-1.45960311115694966380 

-1.45960311115694966350 

-3.000E(-19) 

7.771561(-15) 

1.00 

-1.71828182845904523540 

-1.71828182845904523500 

-4.000E(-19) 

1.065814(-14) 


Table 2 
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Figure 1: Absolute Stability Region 


Using Mat Lab software, the absolute stability region of the new 
method is plotted and shown in figure 1. 

Numerical Example 

Problem I: We consider a highly stiff problem 

/'+1001/+1000;;, j(0) = l,/(0) = -1 


Exact Solution: = exp(-x) h = — 

Problem II: 

f(x,y,y') = y', j(o) = i, /(0) = ^ ,o<x<i. 

Exact Solution: v ( xl = 1 — e* with h =- 

• V ’ 100 

Conclusion 

It is evident from the above tables that our proposed methods 
are indeed accurate, and can handle stiff equations. Also in terms of 
stability analysis, the method is A-stable 

Comparing the new method with the existing method, the result 
presented in the tables 1 and 2 shows that the new method performs 
better than the existing method, and even the order of new method 
is higher than the order of the existing method [11,18], In this article, 
a one-step block method with three off-step points is derived via the 
interpolation and collocation approach. The developed method is 
consistent, A-stable, convergent, with a region of absolute stability 
and order Ten. 
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